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Abstract — A new type of states in graphene-based planar heterojunctions has been studied in the envelope 
wave function approximation. The condition for the formation of these states is the intersection between 
the dispersion curves of graphene and its gap modification. This type of states can also occur in smooth 
graphene-based heterojunctions. 
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Graphene is a promising material for future carbon 
nanoelectronics. Owing to its unique electronic prop- 
erties, this material has attracted particular attention 
of both theoreticians and experimentalists. For exam- 
ple, the mobility of charge carries in graphene can be as 
high as 2-10 5 cm 2 /Vxs and the transport in submicron 
samples can be ballistic [1, 2]. 

We consider a planar hetero junction composed of 
graphene and a gap modification of graphene. When 
we say a gap modification of graphene we imply a 
graphene with an energy gap in the Dirac spectrum 
of charge carriers. There are several gap modifications 
of graphene. 

First, the energy gap can open because graphene 
sheets are located not on Si02 substrate but on some 
other material, for example, hexagonal boron nitride 
(h-BN), when two triangular sublattices of graphene 
become nonequivalent and a gap modification of graphene 
is formed with an energy gap of 53 meV [3]. Second, 
the energy gap opens in the epitaxially grown graphene 
on the SiC substrate [4], which is equal to 0.26 eV ac- 
cording to experimental results obtained by angular- 
resolved photoemission spectroscopy [5]. Third, re- 
cently another modification of graphene, i.e., graphane, 
was synthesized by hydrogenation [6], which has a di- 
rect energy gap of 5.4 eV at the Y point according to 
the calculations [7] . In the first two cases, a graphene 
film deposited on inhomogeneous Si02-h-BN or Si02~ 
SiC substrates can be used (Fig. la shows the case 
with h-BN). In the third case, an inhomogeneously hy- 
drogenated graphene is used (a part of the graphene 
sample is left without hydrogenation, Fig. lb). 

We assume that the energy gap in the gap modifica- 
tions of graphene opens at K and K' points of the first 
Brillouin zone, which correspond to the Dirac points of 
gapless graphene. 

Let us assume that the x axis is directed along 
the plane of the heteroj unction perpendicular to the 
boundary between graphene and its gap modification 



and the y axis is directed along the boundary. The z 
axis is directed perpendicular to the plane of the het- 
eroj unction. The half-plane x < is occupied by the 
gap-less graphene and the half-plane x > belongs to 
the gap modification of graphene. 

The equation for the envelope wave function that 
describes charge carriers in the planar graphene-based 
heteroj unction is written in the form [8, 9] 



[VFj (To ® 0~ x p x 
+ To ® O~o (Vj ~ 



VT Z ®CTyPy) 

E)] (a, y) -- 



- To 



) a z Aj 



(1) 



Here, the parameters with j = 1 are related to the 
gapless graphene and the parameters with j = 2 are 




Fig. 1. Two variants of the system in question: 
(a) graphene layer on the substrate composed of 
h-BN and Si02 and (b) nonuniformly hydrogenated 
graphene on the Si02 substrate. Open circles are 
hydrogen atoms, which are located so that they are 
bound to carbon atoms of one sublattice on one side 
of graphene sheet and to carbon atoms of the other 
sublattice on the other side. 
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Fig. 2. Graphene hetero junction under consideration. 

related to the gap modification of the graphene: vpi 
and Vp2 are the Fermi velocities (in the general case, 
vf2 7^ vpi: and vpi ~ 10 8 cm/s); Ai = and A2 ^ 
are the half- widths of the energy gaps; V\ and V2 are 
the work functions (V2 determines position of the mid- 
dle of the energy gap for the gap modification of the 
graphene with respect to the Dirac points of the gap- 
less graphene, and V\ = is chosen for the origin, see 
Fig. 2). The Pauli matrices a xi a yi a z and 2x2 unit 
matrix <jq operate in the sublattice space (A and B 
sublattices of graphene hexagonal lattice). The Pauli 
matrices r x , r y , r z and 2x2 unit matrix tq operate in 
the valley space (K and K' points of the first Brillouin 
zone). Sign <g> denotes the Kronecker matrix product. 
The operators p x = —i-^ and p y — are the mo- 

mentum operators (h = 1). 

In order to avoid spontaneous generation of elec- 
tron-hole pairs, we assume that the hetero junction in 
question is a junction of the first kind, i.e., the Dirac 
points of gapless graphene are located inside the en- 
ergy gap of its gap modification. This limits value of 
the work function |T^| < A2. 

Motion of charge carriers along y axis is free: 



^(x,y) = ^(x)e ikyV . 



(2) 



The wave function ^f(x) is a bispinor 



where ^k(x) and ^k'{x) spinors describe charge car- 
riers in the K and K' valleys, respectively: 



V K (x) = 



^ka(x) 
^kb(x) 



^K'(x) = 



^k>a(x) 
^k'b(x) 



Let us consider the parity operator 

P = r z <gX7 , 



(3) 



which is a product of the inversion operator 274 = ir z <g> 
<j z and the operator of rotation by the angle tt about 
the z axis A z = —iro (g) a z . Evidently, the operator (3) 
commutes with the Hamiltonian in Eq. (1). 

Equation (1) is solved within the class of wave eigen- 
functions ^\(x) of the parity operator (3) 







*-i(ar) 



A (4) 



Equation (1) can be easily represented as two 2x2 
matrix equations 

E x t> XK (x), 
= E x ^xk'{x). 

(6) 

In this case, we have A = +1 in Eq. (5) and A = -1 in 
Eq. (6). 

Clearly, for Aj = and Vj = 0, we return to the 
spinor wave functions that describe the chiral states ei- 
ther ar the K point or at the K' point. In this case, it is 
possible to introduce helicity operator h = cr • p/ (2 |p|). 
Its eigenvalue (helicity) determines the attribution of 
charge carriers to one of two valleys [10]. However, for 
A / 0, the chiral symmetry is broken, and, therefore, 
instead of the helicity the quantum number A (parity) 
is introduced, which, according to Eq. (4), determines 
attribution of charge carriers to one of the two valleys. 

We use the following condition of matching the en- 
velope wave functions [11, 12] 



(7) 



where the signs "— " and indicate the quantities 
related to the material on the left-hand and right-hand 
sides of the interface, respectively. 

The solution to Eq. (5) for boundary states has the 
form 



^xk(x) 



C(l) exp(^x), x < 0, 
C (qb) exp(-^ 2 ^), x > 0, 



(8) 



where 



.VFljky ~ ttl) 

1 E X " 



. VF2(k y + n 2 ) 
l E x -V 2 + \A 2 ' 



is the constant 

V F 2 



C is the normalization factor, b - 

obtained when matching solutions for x < and x > 
at the line x = under condition (7), 



E x = ±v F1 J kl 



(9) 
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from which it follows that the necessary condition for 
the existence of the boundary states is given by in- 
equality 



«i < \k y \- 
Equation (9) can be rewritten as 



(10) 



«1 = \l tf, 



XI U F1> 



Therefore, the following inequality should also be valid 



< v F1 \k y \. 

Expression for n 2 is represented in the form 
1 



(11) 
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Moreover, the matching leads to the inequality 

VFi(k y - Ki) _ v F2 (k y + n 2 ) 



E x -V 2 + AA 2 ' 



The solution to Eq. (6) is produced from Eq. (8) by the 
following substitutions in factors a and q: k y — ^ —k y 
and A — >> —A. 

Let us discuss separately the case of zero mode 
E\ = 0. Components of the envelope wave function in 
x < region (gapless graphene) ^xk — (^) exp(^ix) 
satisfy equations: 



k y )ai 
k y )a 2 



i.e., either k>i = k y (k y > 0) and a 2 = 0, or k\ = —k y 
(k y < 0) and a\ = 0. Then it follows from the match- 
ing condition (7) that both components of the envelope 
wave function are zero in x > region (b = 0); there- 
fore, we have a\ — and a 2 = 0, i.e., ^xk(x) = 0. 
Thus, there is no zero mode for the boundary states in 
question. 

The following equations are easily obtained from 
Eq. (12): 

E X (E X -V 2 ) _ k2 (13) 



XA 2 E X 



v F1 v F2 
v F1 v F2 k y (^ 1 



(14) 



The two latter equations are valid for either value of A 
(for both valleys), because they are invariant in respect 



to simultaneous substitutions k y — » 



-k y and A — » —A. 



Since k>i > and n 2 > 0, right-hand side of Eq. (13) 
should be positive. Let us denote by So(k y ) such value 
of E\ that the right-hand side of Eq. (13) turns zero, 



holds for hole boundary states (E\ < 0). The bound- 
ary states are not degenerate in parity. That means 
that there is no Kramers degeneracy of energy spec- 
trum for them. This is also true for boundary states 
in a planar quantum well based on graphene nanorib- 
bon [13] and for boundary states localized on zigzag 
edges of gapless graphene [14]. Since parity determines 
charge carrier attribution to one of two valleys, the 
property mentioned above means also that there is a 
"valley polarization" of boundary states: electrons that 
move along the heteroj unction boundary with k y > 
are located near K point and electrons with k y < 
are near K' point and vise versa in case of holes. Be- 
cause of that, current that flows along the heteroj unc- 
tion boundary would be "valley-polarized" . 

By squaring Eq. (14) we get a quadratic equation, 
solution of which produces dependence of energy on k y : 



(12) E X (ky) 



v F iv F -ky]V 2 



Xv F ik y A 2 



h2 



V 2 



Z,2 



(17) 

where v F - = v F i — v F2 . Equation (17) takes into ac- 
count that sign of \k y determines type of charge carri- 
ers in the boundary states. 

It is easy to verify that inequality (11) is always 
true if the energy is given by Eq. (17). Therefore, 
inequality (10) also holds. 

Now, it is simple to analyze inequality (16). Let us 
introduce the following notation: 



k y i — 



\VF- 



lv F2 V% + 2v F _A| =p \ V 2 \yJv F2 V% +4v F1 v F -&l 



2v F2 v 2 F _ 



Under the condition 



vfi < v F2 < 2v FU ^/v F ^\v F ~^A 2 < \V 2 \ < A 2 , 

v F2 

the boundary states exist in the ranges 1 

< \ky\ < ky 2 , kyS < | ky | < kyl 

either for electrons, if V 2 < 0, or for holes, if V 2 > 0. 
Under condition 



vfi < v F2 < 2v FU \V 2 \ < y/v^vp-^ 

VF2 

the boundary states exist in the range 
< \k y \ < k yl 



£o(k y ) 



V 2 



± 



V2 2 



v F1 v F2 k 2 , 



either for electrons, if V 2 < 0, or for holes, if V 2 > 0. 

(15) 

Under the condition 



where "+" corresponds to electrons and "-" to holes. 
Then, the condition ^1^2 > is equivalent to the in- 
equality 

\Ex\ > \eo(k y )\. (16) 



vfi > v F2 , < V 2 < A 2 
the electron boundary states exist in the range 
ky3 ^ I ky I <C kyi , 



It follows from Eq. (14) that inequality Xk y > i Her e and below, we exclude the point k y = 0, because it 
holds for electron boundary states (E\ > 0), and Xk y < corresponds to E x = 0. 
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and the hole boundary states exist in the range 

%2- 



(a) 



< \ky\ < ky 



Under condition 

v F i > v F2 , -A 2 < V 2 < 
the electron boundary states exist in the range 

< \ky\ < ky2, 

and the hole boundary states exist in the range 

ky3 < \ky\ < kyl. 

Let us consider three special cases. 

(1) Under condition V 2 = and vp- ^ 0, the bound- 
ary states exist for both electrons and holes in the fol- 
lowing range if vpi > vp 2 

o<|fej< 



(2) Under condition vfi = VF2, < | | < A 2 the 
boundary states exist in the range 



o< \U < 



VF2\V 2 \ 



either for electrons, if V 2 < 0, or for holes, if V 2 > 0. 

(3) Under condition vpi — VF2-, V 2 = 0, the bound- 
ary states are absent both for electrons and holes, be- 
cause \E\(k y )\ = \so(k y )\, which is in contradiction 
with inequality (16). 

Fig. 3 shows dispersion curves E^ h (k y ) and 6Q h (k y ) 
for the electron and hole boundary states for three val- 
ues of V 2 in the model of graphene-based heteroj unc- 
tion with A 2 = 260 meV and vp 2 = 1.2 x 10 8 cm/s 
for gap modification of graphene. 

Our results remain in essence the same if instead of 
a sharp heteroj unction we consider a smooth het ex- 
junction. Indeed, let vp(x) and A(x) vary smoothly 
from their values for gapless graphene to their values 
in gap modification of graphene over a strip with the 
width d < k>i 1 2 . Then change in energy of the bound- 
ary states is \SE\(k y )\ < 1 meV. Such insignificant 
variation in energy of the boundary states produces 
no noticeable qualitative changes. A similar result has 
been obtained for boundary states in heteroj unctions 
of narrow-gap semiconductors with intercrossing dis- 
persion curves in [15]. 

To conclude, we would like to point out that the 
new type of boundary states in graphene heteroj unc- 
tions can be studied in experiment by tunnel spec- 
troscopy of angular-resolved photoemission spectroscopy 
similar to how it have been done for boundary states 
in gapless graphene [16-18]. 
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(a) there are no boundary states for electrons and 
holes at V 2 = 0, (b) there are only hole boundary 
states in the range < \k y \ < k y \ at V 2 — 100 meV, 
and (c) there are only hole boundary states in the 
ranges < \k y \ < k y2 and k y s < \k y \ < k y \ at V 2 = 
250 meV. 
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